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ABSTRACT 

We report a preliminary result for D°-D° mixing and the doubly Cabibbo 
suppressed decay rate Rd based on an analysis of D° — * K + ir~ decays from 
57.1 fb _1 of data collected at or just below the T(AS) resonance with the 
BABAR detector at the PEP-II collider. We set 95% confidence limits for the 
mixing parameters x and y' and find that our result is compatible with no 
mixing and no CP violation. In the limit of no mixing we find the doubly 
Cabibbo suppressed decay rate R D = (0.357 ± 0.022 (stat.) ± 0.027 (syst.))% 
and the CP violating asymmetry A D = 0.095 ± 0.061 (stat.) ± 0.083 (syst.). 



1 Motivation 



Mixing can be characterised by the two parameters x = Am/T and y = Ar/2r, 
where Am (AT) is the difference in mass (width) between the two different mass 
eigenstates and T is the average width. 

* On behalf of the BABAR collaboration. 



Within the Standard Model the level of D°-D° mixing and CP violation 
is predicted to be below the sensitivity of current experiments 0) . For this 
reason it is a good place to look for signals of physics beyond the Standard 
Model. Other experiments ID have already tried this with smaller datasets 
using a technique similar to what is described here. In any attempt to measure 
mixing one should consider the possibility of CP violation also as, with new 
physics, there is no a priori expectation that it is insignificant. 

Mixing and CP violation can be detected by observation of the wrong-sign 
decay D° — > K + ir~ (charge conjugation is implied unless otherwise stated). 
Production through direct decay is doubly Cabibbo suppressed (DCS) but it is 
also possible for the D° to oscillate into a D° and subsequently decay through 
the right-sign Cabibbo favoured (CF) decay D° — > K + tt~ . The two processes 
can only be distinguished by an analysis of the time evolution of the decay. 

Assuming a;', j/' C 1 and that CP is conserved, the time-dependent decay 
rate for the wrong-sign decay D° — > K + ir~ from DCS decays and mixing is 
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where r^o is the D° lifetime and Rd is the ratio of DCS to CF decaysQ Because 
x' only appears in the time distribution as a squared value, it is not possible to 
determine the sign of x' in an analysis based on the D° — > K + ir~ decay alone. 

CP violation can be either direct, in mixing or in the interference between 
the two. The CP violation gives rise to different apparent values for the pa- 
rameters in eq. |l| so we define -Rp^~\ x' + ^ and y' + ^ for D mesons produced 
as a D° (D°). 

2 The BABAR detector and data selection 

For this analysis, we use 57.1 fb _1 of data collected with the the BABAR de- 
tector which is described in detail elsewhere P. Reconstruction of charged 
particles and particle identification of Kaons and pions are the most essen- 
tial. Tracking is provided by a five-layer silicon vertex tracker (SVT) and a 

1 x' — x cos Sk-k + U shi 5ktt and y' = —x sin Sk-k + V cos 5kk where Sk-k is an 
unknown strong phase. 



forty-layer drift chamber (DCH), both in a 1.5 T solenoidal magnetic field. A 
Cherenkov ring imaging detector (DIRC) is placed outside the tracking volume. 

Kaons (pions) are identified by calculating a likelihood product of the 
information from dE/dx measurements in the SVT and DCH and the recon- 
structed Cherenkov angle and photon statistics from the DIRC, with an effi- 
ciency above 75% (80%) and mis-id rate below 8% (7%) for p < 4 GeV/c. 

We select D° candidates from reconstructed D* + — > D°tt + decays. The 
charge of the pion in the decay identifies the flavour of the D° and also serves 
to create a clean sample of D° decays. Both right-sign and wrong-sign D° 
candidates are selected. We select only D* + candidates with Pd-+ > 2.6 GeV/c 
in the centre-of-mass frame to reject D* + candidates from B decays. Other 
event selection criteria are employed to ensure that we have high quality tracks 
and do not have any D* + candidates with multiple overlapping tracks. 

3 Results 

An unbinned maximum likelihood fit is used to extract the mixing parame- 
ters. For each D* + candidate we use the D° candidate mass mjf„ the mass 
difference Sm between the D* + and the D° candidate and the proper lifetime 
and error on the lifetime of the D° candidate. The D° lifetime and the signal 
resolution model is determined from the large right-sign sample. Sidebands are 
included in vhk-k and Sm such that the level and time evolution of the different 
background types can be evaluated. In figure [I] we show the time evolution 
of the wrong-sign sample. In total we observe around 120,000 (440) right-sign 
(wrong-sign) signal events. 

Since the fit allows x' 2 to take unphysical negative values an error estimate 
from the log-likelihood surface (LLS) would require a Bayesian analysis where 
the choice of prior is not clear. In addition, an accurate error estimate from 
the LLS requires a LLS shape that is independent of the outcome of the fit. 
At the current level of statistics these requirements are not even approximately 
met, especially for the small mixing values observed. 

Instead we use a method where we define a 95% confidence limit contour in 
x' 2 and y' space using toy Monte Carlo experiments^]. Contours are constructed 

2 With a toy Monte Carlo experiment we mean a Monte Carlo sample of the 
same size as the data generated from the PDF of the fit. 



such that there is a 95% probability for any point a c = (x' c , y' c ) on the contour 
that the likelihood ratio 

Aln£(a c ) = ln£(<3 c ) - ln£ max , (2) 

will be greater than the corresponding value Aln£data(o?c) calculated for the 
data. £ m ax is here the maximum likelihood obtained from the fit to either data 
or a toy Monte Carlo sample. The probability is evaluated by creating multiple 
toy Monte Carlo samples at the point a c and for each of the samples evaluate 
Aln£(d< c ) after a fit. 

As well as for the general case allowing for CP violation we also calculate 
our results for the special cases where CP is conserved and where no mixing 
is allowed. In the case where we assume no mixing we calculate the direct CP 
violation term Ad = (R^ — R^)/(R^ + R^). 

Table 1: A summary of our results including systematic errors. A central value 
is reported for both the full fit allowing x' 2 < 0, and from a fit with x' 2 fixed at 
zero. The 95% confidence limits are for the case where x' 2 is free. 



Fitted Central Value 



Fit type 


Parameter 


x' 2 free 


x' 2 fixed at 


95% C.L. interval 




K [%] 


0.32 


0.35 


0.18 <R&< 0.62 






0.26 


0.27 


0.12 <R D < 0.56 


CP violation 


x'+ 2 


-0.0008 





x'+ 2 < 0.0035 


allowed 


x '-2 


-0.0002 





x'- 2 < 0.0036 




y' + [%] 


1.7 


0.7 


-7.5 < y'+ < 3.4 




y'- [%] 


1.2 


0.9 


-5.7 < y'~ < 3.6 


No CP 


Rd [%] 


0.30 


0.31 


0.22 < i? D < 0.46 


violation 


x' 2 


-0.0003 





x' 2 < 0.0021 




y' [%} 


1.3 


0.8 


-3.7 < y' < 2.4 



No mixing R D = (0.357 ± 0.022 (stat.) ± 0.027 (syst.))% 
A D = 0.095 ± 0.061 (stat.) ± 0.083 (syst.) 



The confidence contours for the mixing results including systematic errors 
are shown in figure |^ and the overall results are summarised in table [l| 

For our systematic errors we evaluate the contributions from uncertainties 
in the parametrisation of the PDF's, detector effects, and effects of our selection 
criteria. For detector effects like alignment errors or charge asymmetry we 
measure their effect on the right-sign control sample. For variations in the 
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Figure 1: A projection of the wrong- Figure 2: 95% confidence limits for 

sign data on the time axis with the D° (dashed) and D° (dotted) sepa- 

main plot showing the signal region rately (allowing for CP violation) and 

and the inset a sideband. The plots the contour when no CP violation is 

show data as points and the resulting allowed (full). The solid point repre- 

projection of the background (shaded) sents the most likely fit point in the 

and signal (open) from the fit. case of no CP violation and the open 

circle the same but forcing x' 2 > 0. 

event selection we assign for this preliminary result the full variations in the 
resulting contours as systematic errors. 

In summary we have set new and improved limits on mixing and CP 
violation for neutral D mesons. Our results are compatible with no mixing and 
no CP violation, all of which fits well with the predictions from the Standard 
Model given our current sensitivity. 
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